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§1. Three-parameter Free-damped Modes 

In 1998, Rosu and Reyes 1 ) introduced classical three-parameter free damped 
oscillators by means of the nonrelativistic supersymmetry (factorization) technique. 
They showed that in this approach the usual two-parameter free damped oscillator 
equation (the dots are for time derivatives) 

y" + 2py + u%y = 
is generalized to a damped oscillator equation of the form 

y + 2(3y + uj\t; i)y = , where u 2 (t; 7) = Lg - ^ T^ l ' ( M ) 

possessing damping solutions with singularities that we call here 7 modes, see below. 
The new damping solutions depend on an additional parameter 7 given by the general 
Riccati solution associated to the original two-parameter damped equation. This 7 
parameter introduces a new time scale in the damping phenomena that is related to 
the time at which the singularities occur. 

For the three types of free damping, the 7 modes are the following: 



(i) Underdamping: j3 2 < uj u = \/cj'q — ft 
Vu — ~A u e ^ 



7 

cj u sm(uj u t + (/)) + — — — - cos(cj u t + (/)) . (1-2) 
7c + 1 



typeset using pTpTEX.ds (Ver.0.9) 
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(ii) Overdamping, (3 2 > Uq, lj = \f f3 2 — cJq, 



Vo -A Q e 
(iii) Critical damping, fS l — cJq. 



cj g sinh(cj t + 0) cosh(cj £ + 0) . (1-3) 

jt + 1 



7* + 1 7 2 



+ -|(7t + l) 2 • (1-4) 



Damped motion of this type could be quite common in nature, although if the 
value of 7 is very small it is difficult to distinguish it from the usual two-parameter 
free damping. In such cases it is still possible to estimate 7 through some powerful 
estimation methods. The problem of parameter estimation in dynamical systems has 
motivated extensive studies of the adaptive observer designs for linear and nonlinear 
systems during the last decade. Among the single output adaptive schemes that can 
be directly applied to our case we cite those of Marino, 2 ) Marino and Tomei, 3 ) and 
Besancon. 4 ) It is the purpose of this paper to show that very good estimates for 7 
can be obtained by applying the latter scheme. The main reason to use the Besancon 
scheme is its simplicity with respect to the scheme of Marino that requires solving 
more differential equations in real time in order to obtain the corresponding velocity 
to recover the unknown value of the parameter (s). In the next section, Eq. (|l-ip is 
written as a system of two first-order differential equations in the afflne form which 
is required in order to apply Besancon's adaptive observer scheme. This is detailed 
in Section 3 where we also give some pedestrian definitions of the technical concepts 
involved in the adaptive estimation algorithms. A small conclusion section ends up 
the paper. 

§2. Equivalent System of First Order Differential Equations 

Denoting X\ — y and X2 = y, Eq. (|l-ip can be rewritten as 

Xi = X 2 (2-1) 
X 2 = -2(3X 2 + {j^^ ~ 4) Xi . (2-2) 

We will consider X\ as the naturally measured state (the most easy to measure). 
Therefore, it seems logical to take X\ as the output of the system, that is y = 
CX, where C = [1 0] and X is state vector. In addition, we consider that (3 is 
available in some way and only 7 stands as an unknown parameter. Thus, the 
complete monitoring of the system means the on-line estimation of X2 and the on- 
line identification of the 7 value. The available methods for parameter identification 
need that the unknown parameters are in afflne form, i.e. 

Xi = X 2 (2-3) 
X 2 = -2/3X 2 + aXx , (2-4) 
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where a is the unknown constant parameter. It is clear that 7 is not an affine 
parameter but we can try to estimate the time-dependent quantity 



because it appears in the affine form. Notice that, even though a is not a constant, 
its evolution is significant only close to the origin, that is, as time increases a turns 
rapidly into a constant parameter. 



It is well known that only in very special occasions one can have a sensor on 
every state variable, and some form of reconstruction from the available measurement 
output data is needed, in many cases. In general, an observer is an algorithm that 
reconstructs the unmeasurable states of a system from the measurable output. The 
so-called high gain techniques proved to be very efficient for state estimation, leading 
to the well-known concept of high gain observer.^ The gain is the amount of increase 
in error in the observer's structure. This amount is directly related to the velocity 
with which the observer recovers the unknown signal. The high-gain observer is an 
algorithm in which the amount of increase in error is constant and usually of high 
values in order to achieve a fast recover of the unmeasurable states. However, when 
a system depends on some unknown parameters, the design of the observer has to be 
modified in such a way that the state variables and parameters could be estimated. 
This leads to the so called adaptive observers, i.e., observers that can change in order 
to work better or provide more fit for a particular purpose. 

Here, we present an interesting application of the recent scheme of Besancon and 
collaborators, 4 ^ for unknown parameters that appear linearly in the dynamical sys- 
tem to the interesting case of state and parameter estimation of the three-parameter 
free-damped oscillator. In this particular case, the unknown parameter appears in 
a nonlinear manner and in order to avoid this difficulty we group together a set of 
parameters which leads to a time-dependent parameter. The block scheme of the 
adaptive observer for this case is presented in Fig. (1). 

The main assumptions on the considered class of systems are that if all param- 
eters were known, some high-gain observer could be designed in a classical way, and 
that the systems are "sufficiently excited" in a sense which is close to the usually 
required assumption on adaptive systems (the signals must be rich enough so that 
the unknown parameters can indeed be identified). 

We write the 7 free damping motion system given by Eqs. (|2-3I) - (I2-4I) in Be- 
sancon's form: 



where: 

A - Brunovsky-type matrix, i.e., of the form A = where 5ij is the Kronecker 

symbol 




(2-5) 



§3. Design of the Adaptive Observer 



F = 



x = ax + tp(x) +&(x)e 

y = CX 
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X n GK - n-component real vector 
y - measured output 

\P(t) E W nXp - the vector of known functions (p = the number of unknown parame- 
ters) 

9 E P - vector of unknown parameters (i.e., a vector belonging to some known 
compact set that should be estimated through the measurements of the output y). 

For the 7 system, we have: n — 2, p = 1 and 9 — a. 

Therefore the entries of the matrix A are an = 0, ayi = 1, &21 = and (22,2 = 0, the 
entries of the vector \P are ^11 = and ^21 = whereas the entries of the vector 
(f are <pn — and <p2i — —2/3X2- 



We consider now the following dynamical system 



' j£ = AX + + V(a(X))a(6) + [pTf + rr T C T ] (y - CX) 

§ = [ p ^r T c T ] (y -ex) 

t = p{A-KC)r + pV(a{X)) , 




where r E M nXp , K is some vector for which A — is a stable matrix, A — 
diag[l, p _1 , . . . , p~( n_1 )], p E R+ being a constant to be chosen, and a is a so-called 
saturation function. The saturation function is a map whose image is bounded by 
the chosen upper and lower limits, B and 6, respectively. 6 ) This avoids the over 
and/or under estimation, keeping the possible values within the set of the physically 
reliable values, and consequently increasing the chance of the quick convergence to 
the true values. Thus, the latter is the following map 



a(s) = 



F has been proven by Besancon to be a global exponential adaptive observer for 
the system F . That means that for any initial conditions X(£q), X(t$), 0(to) and 
G F, the errors X(t) — X(t) and 6(t) — 9{t) tend to zero exponentially fast when 
t —> 00. Consequently, for the free damping motion system, the matrix A — KC has 
the eigenvalues 

\ 12 = -1/2 fei ± 1/2 ^h 2 -4k 2 (3-1) 

Selecting k 2 = 5 we the eigenvalues Ai = A2 = — fci/2, and choosing any 

k\ > we make A — KC a stable matrix. Thus, the explicit non matrix form of the 
observer given by F is the following 



x 1 = x 2 + ( P h + a 2 ) (x x - x{) 
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X 2 = -2pa(X 2 ) + <t(Xi)<7(0) + p + AAJ (X 1 - Xi) 

A = p(-fciA + A) 

A = -^i 2 A + p<7(X 1 )) . 

To recover now the value of the parameter 7 from the estimated value of a, we solve 
for 7 the equation 2(7/(7^ + l)] 2 — Uq = a, giving 



-t(at + /3 2 )- ^2a + 2^ 
7 = 1?(a + p)-2 • (3 ' 2) 

We still can hope that the value of 7 can be estimated with sufficiently good accuracy 

until the blow-up time td = y/2/ [a + /3 2 ) 2 and indeed this is what we have found 
for the 7 modes. 

In the computer simulations that we performed for illustrating purpose we took 
always p — 100 for the constant parameter defining the diagonal matrix A. Fig- 
ures (2) - (5) show the results of numerical simulations for an underdamped case 
with the following parameters: (3 — 0.04, 7 = 0.15, k — 1.00, uo = 1.00. In figure 
(6), we present numerical simulations for the Riccati parameter 7 for the overdamped 
case with (3 = 1.0 and uo = 0.2 in the plot (a) and the critical case with (3 = uo = 0.1 
in the plot (b). Figure 5a shows the initial stages of the process of convergence dis- 
playing big oscillations. After this, one can always notice the convergence to the 
true value once the estimate enters the neighborhood of the true value determined 
by the exponential criterium of convergence to the neighborhood, where small os- 
cillations around the true value are quite natural. The numerical simulations have 
been performed with the relative values of the parameters chosen to illustrate the 
three classes of damping, but otherwise there are no restrictions on their values. 



§4. Conclusion 



The adaptive observer scheme of Besancon et al is an effective procedure to 
rebuild the Riccati parameter of the class of damping modes of Rosu and Reyes. De- 
spite the fact that in the estimation algorithm we used an associated time-varying 
quantity for which the exponential convergence is not guaranteed, we obtain excellent 
results for the estimates of the unknown constant Riccati parameter. As it can be 
seen from this work, the application of Besancon's adaptive scheme is not necessarily 
restricted to the estimation of constant parameters but one can go with it to some 
time-varying cases. In addition, the scheme provides very good estimates for the 
unknown states of this class of oscillators for each of the three possible cases. More- 
over, the rapid identification of the additional parameter 7 of the chirped damped 
oscillators by this method opens a way to the experimental study of these oscillators. 
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Fig. 1. Schematic representation of the adaptive observer, where the output of the system is the 
input of the adaptive observer and the outputs of the latter are the rebuilt states and the 
parameter a from which it is possible to rebuild 7. 




time time 
Fig. 2. Numerical simulations for the state X± in the underdamping case. In (a) the solid line 
represents the time evolution of the true state and the dotted line gives the evolution of its 
estimate. The plot (b) is a detail of the figure (a) to appreciate the variation of X\ in the 
beginning. 




time time 
Fig. 3. Numerical simulation for the state X2 in the underdamping case. In (a) the solid line 
represents the time evolution of the true state and the dotted line gives the evolution of its 
estimate. Plot (b) is a detail of (a) to see better the variation of X2 in the beginning. 
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Fig. 4. Numerical simulation for a in the underdamping case. The solid line represents the true 
value of a and the dotted line gives the estimate. Plot (b) is a detail of (a) to see better the 
variation of the estimated a in the beginning. 
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Fig. 5. Numerical simulation for the parameter 7 in the underdamping case. The true value of 7 

is given by the solid line whereas the dotted line gives the estimate. Plot (b) gives the details 
of (a) near the origin to better see the variation of the estimate of 7. 
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Fig. 6. Numerical simulations for the parameter 7 in the overdamped case [plot (a)], and in the 
critical case [plot (b)]. The solid line gives the true value of 7 and the dotted line gives its 
estimate. 



